Saffman- Taylor streamers: mutual finger interaction in an electric breakdown 
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Bunches of streamers form the early stages of sparks and lightning but theory presently concen- 
trates on single streamers or on coarse approximations of whole breakdown trees. Here a periodic 
array of interacting streamer discharges in a strong homogeneous electric field is studied in PDE 
approximation in two dimensions. If the period of the streamer array is small enough, the streamers 
do not branch, but approach uniform translation. When the streamers are close to the branching 
regime, the enhanced field at the tip of the streamer is close to 2_E tx) , where Eoo is the homogeneous 
field applied between the electrodes. We discuss a moving boundary approximation to the set of 
PDEs. This moving boundary model turns out to be essentially the same as the one for two-fluid 
Hele-Shaw flows. In two dimensions, this model possesses a known analytical solution. The shape of 
the 2D interacting streamers in uniform motion obtained from the PDE simulations is actually well 
fitted by the analytically known "selected Saffman- Taylor finger" . This finding helps to understand 
streamer interactions and raises new questions on the general theory of finger selection in moving 
boundary problems. 



I. INTRODUCTION 

Streamers are growing ionized fingers that appear in 
electric breakdown whenever non-ionized matter is sud- 
denly exposed to strong electric fields, therefore they are 
very common in nature and technology in gases, liq- 
uids and solids [H, H, [1]. They occur for instance in 
early stages of atmospheric discharges such as sparks 
and lightning or in sprite discharges high above thun- 
derclouds 0, IE IE S m ■ Streamers are characterized by 
a thin space charge layer around their tip that enhances 
the local electric field; this enhanced field in turn creates 
a very active impact ionization region. 

Most experiments produce many streamers, certainly 
when the emitting electrode is a long wire [9|, a nd fre- 
quently also when it is the point of a needle [10| . Sim- 
ulations, on the other hand, concentrate almost exclu- 
sively either on single streamers within a microscopic dis- 
charge model, or on the complete streamer branching tree 
in quite coarse phenomenological models. Only in pd| . 
the electrostatic interaction of narrow streamers within a 
widely spaced streamer array is studied in relatively low 
electric fields within a microscopic model; as the streamer 
radius is fixed, the numerical implementation is essen- 
tially one-dimensional. In the present paper, we mimic a 
similar periodic array of identical parallel streamers, but 
in a higher field, see Fig. [TJ Furthermore, rather than 
fixing radius and shape of the streamers a priori, we let 
it emerge dynamically within the simulation. Such ar- 
rays of streamers can be created experimentall y b y an 
array of needles inserted into a plate electrode [12l [l3| . 
Bunches of parallel streamers have also been observed in 
sprite discharges above thunderclouds [1, Q . 

The problems addressed in this paper are multiple: 
What is the charge distribution and velocity of an ar- 
ray of streamers, depending on their distance and on the 
applied electric field? Do they approach a state of uni- 
form translation, in contrast to single streamers? And 



FIG. 1: Periodic array of negative streamers (net charge den- 
sity) in a strong homogeneous background electric field Eoo 
pointing downwards. L is the period of the array. The dash 
lines represent two symmetry lines. The box around a part of 
the central streamer indicates the part presented in Fig. [3j 



how can the dynamical evolution of their shape be placed 
in the context of other moving boundary problems in 
nature? Giving already a major conclusion of the pa- 
per, we find that uniformly translating streamer arrays 
in the microscopic discharge model in two spatial dimen- 
sions are very well fitted by a classical solution of two- 
fluid-flow na mely by the so-called selected Saffman- 
Taylor finger [15| , cf . Fig. [31 Therefore the velocity of the 
streamer array is about twice the electron drift motion 
in the background field, and their diameter approaches 
half the period of the array. The observation also raises a 
theoretical question on pattern selection, namely why the 
same finger shape is selected in the hydrodynamic and in 
the discharge problem, given the fact that the problems 
are similar but not identical. 
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The paper is organized as follows: Sections III Al and 
III Bl introduce the minimal PDE model for streamers 
adapted to describe the evolution of an array of stream- 
ers. The general behaviour and properties of these in- 
teracting streamers are discussed in Sec. Ill CI Section 
IIII Al presents a moving boundary approximation of the 
minimal model used in the simulations. In Sec. HUB I we 
present an analytical solution of this approximation for 
the shape of the streamer. This solution is known as the 
selected Saff man- Taylor finger and fits well the charge 
distribution of the streamer, cf. Fig. [3] In Sec. IIII CI we 
briefly discuss some open issues related to this boundary 
observation. We conclude by shortly summarizing and 
discussing our study in Sec. IIVI 



II. 



DENSITY APPROXIMATION AND 
SIMULATION RESULTS 

A. Minimal streamer model 



We analyze negative streamers in simple media like 
pure nitrogen within the minimal streamer model [l], [H], 
all , [HI, [H, H(| that includes electron diffusion and drift 
in a self-consistent electric field, while ions are taken as 
immobile due to their much larger mass. New charge 
carriers are generated by an impact ionization term in 
Townsend approximation that depends nonlinearly on 
the local electric field. In dimensionless units, the model 
is 



dtp 
V 2 <f> 



DV 2 a + V- (crE) + er|E|e~ 1/|E| , (1) 
ame-VW, (2) 
G-p, E = -V<j>, (3) 



where a and p are the electron and ion densities, <f> is the 
electrostatic potential, E the electric field, and I? is a 
diffusion coefficient, taken as D = 0.1 The intrin- 

sic length scale of the model is the mean free path of an 
electron between two ionizing collisions in fields |E| ^> 1, 
for nitrogen at standard temperature and pressure, it is 
2.3 pm; the scale of time is 3 ps and the scale of the 
electric field is ~ 200 kV/cm in this case. A general dis- 
cussion of dimensions can be found, e.g., in [U [H, [2l| . 
There it is argued that the main advantage of working 
with dimensionless quantities is that all basic results are 
immediately generalized to any gas pressure, tempera- 
ture and composition. 

The model is solved numerically on adaptively refined 
comoving grids as described in detail in [22|; the finest 
grid in our simulations was 1/4. 

Notice that the only ionization source in our model 
is impact ionization. We assumed this for the sake of 
simplicity and in order to emphasize the elementary pro- 
cesses that are common between streamers and two-phase 
hydrodynamic systems, as will be discussed below. In ex- 
change for this simplicity we restricted ourselves to neg- 
ative streamers in media where photo-ionization is ab- 



sent or negligible, such as pure nitrogen, argon or GaAs. 
Moreover, in 21( it was shown that under certain condi- 
tions the effects of photo-ionization on the propagation 
of negative streamers are negligible even in ambient air. 



B. Implementing an array of streamers 

Another simplifying assumption is the restriction of 
the problem to two-dimensions. Indeed, we believe that 
the main characteristics of the dynamics of the interact- 
ing streamers in two-dimensions, as described below, will 
be qualitatively the same as those in three-dimensions. 
This is supported by past simulations of single stream- 
ers in two-dimensions [23[ which are qualitatively very 
similar to three-dimensional simulations [21(. Moreover 
the generalization to a three-dimensional geometry is 
not straightforward since the numerical implementation 
of streamers in a 2D periodic lattice is non-trivial due 
mainly to the implementation of the boundary condi- 
tions. Another reason to consider only a two-dimensional 
geometry is that it allows the construction of moving 
boundary approximations with a known analytical solu- 
tion. The agreement between that explicit solution and 
the actual shape of the front is remarkable and detailed 
below. However, two-dimensional streamers are not just 
interesting from an academic point of view, they also oc- 
cur in experiments in thin semiconductor wafers [24| . 

An infinite, periodic array of streamers can be reduced 
to the simulation of a single streamer in a channel with 
Neumann conditions on the lateral boundaries. This is 
done as follows. If the streamer is centered at y = and 
propagates along the x direction, and if the period of 
the streamer array is L, there are two symmetry lines at 
y = ±L/2 where all normal derivatives vanish; therefore 
Neumann conditions for potential and electron density 
d y (j) = 0, dy<7 = 0, at y = ±L/2 can substitute the other 
streamers in the array. 

This array of streamers is now studied in a constant 
electric field E = —EooX far ahead of the streamers; this 
field is imposed as an inhomogeneous Neumann boundary 
condition on <f> at the boundary at x 3> 1 while at x — 
the electrostatic potential is fixed. These conditions 
can be used when planar electrodes are first charged and 
then insulated; for more general electric circuits, they 
also approximate streamers that are much shorter than 
the inter-electrode distance. For the particle densities, 
we used homogeneous Neumann boundary conditions on 
all boundaries. As initial conditions we used in this paper 
an electrically neutral Gaussian seed centered at (x,y) — 
(0, 0) of width 16 and height 1/4.7, except in Fig. [5b, 
where the width is larger but the total number of particles 
is not changed. We verified that the same attractor of the 
dynamics is approached after sufficiently long time when 
lateral position, width and height of the initial seed were 
varied. 
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FIG. 2: Simulated streamers that branch (x) or translate 
uniformly (•) as a function of the period L of the array and 
of the uniform field ahead of it. The line interpolates the 
phase transition. 



C. Branching versus uniform translation 

After a transient evolution, the simulated streamers ei- 
ther reach a state of uniform translation, i.e. they propa- 
gate with constant velocity and unchanged shape, or they 
branch like single streamers [25|, |26| . Two parameters 
control the two regimes for the evolution of the streamers: 
the period of the array, L, and the background electric 
field applied between the electrodes, Eoo- Fig. [2] shows 
a phase diagram spanned by the electric field Eoo and 
the spatial period L; here L = 96 to 616 was explored 
in steps of AL = 40 and E^ — 0.4 to 1.0 in steps of 
AEoo = 0.1. Below the transition line, i.e., for small pe- 
riod L, the proximity of the other streamers suppresses 
branching and the whole streamer array propagates uni- 
formly after some transient stage, while above the line 
the streamers branch eventually. We remark that in gen- 
eral, there can be uniformly translating solutions in the 
part of the phase diagram marked as "branching" ; how- 
ever, the set of initial conditions for which those solutions 
emerge (their basin of attraction) is so small that they 
are not reached from our initial conditions. 

We now analyze in detail the uniformly translating 
streamer array that emerges for sufficiently small E^ 
and/or L (the lower part of the phase diagram, see 
Fig. [5]). After initial transients of duration t = 100 or 
less, these streamer heads reach a constant velocity and 
a constant shape for the rest of the evolution: this is the 
attractor of the dynamics, namely, the solution reached 
after a sufficiently long evolution from a large set of ini- 
tial conditions. Therefore this attractor does not depend 
on the particular choice of the initial seed used. Only 
the transient evolution and its duration can depend on 
this choice. But all these various transient regimes lead 
to the same final uniformly translating state, with the 
same shape and velocity. Fig. [3] shows the space charge 
distribution p — a of the attractor for Eoo — 0.5 and 
L = 256 at the time t — 1800 (long after the transient 
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FIG. 3: (Color online) The thin space charge layer p — a 
around the uniformly translating streamer discharge (density 
color coded) with the Saffman- Taylor finger of width L/2 su- 
perimposed (thick solid line). L is the width of the Hele-Shaw 
cell for the Saffman- Taylor finger or the period of the array 
for streamers, the lateral boundaries then being lines of mir- 
ror symmetry between the streamers. Here L = 256, and 
the electric field far ahead is Eoo =0.5; this corresponds to 
~ 0.059 cm and ~ 100 kV/cm for nitrogen under normal con- 
ditions. Equipotential lines are also plotted (thin solid lines). 



evolution ended) . Figure [4] shows the electric field and 
the net charge profiles from the same simulation. The 
electric field along the streamer axis (y = 0) is presented 
at times 1400, 1600 and 1800 together with the net charge 
density at time 1800. The three profiles of the electric 
field show that the propagation indeed is uniform. The 
thin space charge layer creates a strong field enhance- 
ment immediately ahead of the ionization fronts like in a 
freely propagating streamer. However, behind the space 
charge layer, the electric field profile inside the streamer 
array shows characteristic differences to the field profile 
within a single streamer [HlH^HSHfJ, 26]. Immediately 
behind the space charge layer, the electric field decays 
very rapidly like in a single streamer. Then a transition 
to a slower field decay sets in. Finally, far behind the 
streamer head, the electric field vanishes completely, in 
contrast to the nonvanishing residual field inside a sin- 
gle streamer. These observations require further studies. 
However, one conclusion can already be drawn by apply- 
ing the Poisson equation V • E = p — a to the streamer 
head front as a whole. As the field has a constant value 
—Eoo for ahead of the streamer array and vanishes far be- 
hind the streamer heads, the streamer heads must carry 
an average charge —Eoo per unit area, i.e., each streamer 
head must carry a total charge overshoot of —Eoo • L to 
collectively screen the electric field completely behind the 
array of heads. 




1200 1300 1400 1500 1600 1700 1800 1900 



FIG. 4: Absolute value of the electric field (solid lines) for 
times t = 1400, 1600, and 1800, and space charge density 
(dotted line) at time t = 1800 on the streamer axis for L — 256 
and Eoo = 0.5. Field and density at time t — 1800 correspond 
to the uniformly translating finger in Fig. [3] 



These properties of an array of streamers contrast 
strongly with those of a single streamer, discussed exten- 
sively in [23j . For example, a single streamer in a strong 
homogeneous background electric field never reaches a 
state of uniform translation. The radius of curvature of 
the head of the streamer expands during its motion up to 
the time where instabilities grow and branching occurs. 
Furthermore, the electric field inside of a single streamer 
is not as perfectly screened. There are, therefore, remark- 
able qualitative changes in the propagation of a streamer 
when the interaction of neighbouring streamers is signif- 
icant. These effects are expected to persist also in the 
three-dimensional case. 

However, global considerations on the charge content 
of the streamer head do not fix the shape of the finger 
and the spatial charge distribution within each uniformly 
translating streamer head. These density distributions 
and the consecutive field enhancement and velocity are 
problems of dynamical selection that will be addressed 
in the remainder of the paper. 



III. MOVING BOUNDARY APPROXIMATION 
AND SAFFMAN- TAYLOR SOLUTION 

A. Moving boundary approximation 

As shown in Figs. [1] and [3J after a sufficiently long 
evolution, during the steady evolution of the streamer, 
the width of the ionization front can be much smaller 
than its radius of curvature. Similarly to other pattern 
forming systems, such as solidification fronts, this sepa- 
ration of scales enables one to consider the front as an 
infinitcsimally thin, sharp moving interface. The origi- 
nal nonlinear dynamics is then replaced by a set of linear 
field equations (typically Laplace) on both sides of the 
interface, with appropriate boundary conditions at the 
interface and further away from it. The interface dynam- 
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FIG. 5: Schematic view of the mathematical setup for the 
moving boundary approximation. <f> + and (j>~ stand for the 
electric potential outside and inside the streamer respectively. 



ics is then typically related to gradients of the Laplacian 
fields at its vicinity. 

In 0, 111 S3, El IH, a moving boundary approxi- 
mation was proposed and elaborated for the thin space 
charge layer and associated sharp ionization front that 
separates the ionized from the non-ionized region. This 
model was proposed to describe the evolution of a sin- 
gle streamer but it is straightforward to adapt it to the 
evolution of an array of streamers since we just need to 
add homogeneous Neumann boundary conditions on the 
symmetry lines. See Fig. [5] for a schematic view of the 
mathematical setup of the moving boundary approxima- 
tion. 

The non-ionized, electrically neutral region outside 
the streamer is fully described by V 2 = 0, and <j) — * 
4>o + -E-ooiE fixes the homogeneous field E x far ahead of 
the streamers at x 3> 1. The symmetry line between two 
streamers is represented by a Neumann boundary condi- 
tion for the electric potential, d y 4> = 0, at y = ±L/2. If 
the boundary motion is approximated by the local elec- 
tron drift velocity v = V</>, the interior of the streamer 
as ideally conducting <f> = const, (where the constant can 
be set to due to electrostatic gauge invariance) , and the 
electric potential across the boundary as continuous, we 
arrive precisely at the unregularized moving boundary 
problem for a Saffman- Taylor finger after simply substi- 
tuting the electric potential <fi by the pressure field p. 
This is a classical problem where a very viscous fluid is 
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penetrated by a much less viscous one within the narrow 
spacing of a Hele-Shaw cell. 



B. Comparison with the Saffman-Taylor solution 

An explicit uniformly translating solution for this mov- 
ing boundary problem was found long ago by Saffman 
and Taylor [15[ . The solution for the interface x = x{y, t) 
in a channel of width L is given by 



L(l-A) , 

x = In 

2?r 



1 / /2ttv 
- 1 + cos — - 

2 V V AL 



■vt. 



(4) 



where the velocity is v — £oo/A in our notation and the 
field at the tip is enhanced by a factor 1/A. The param- 
eter A is the ratio between the width of the finger and 
the width of the channel; A can take any value between 
and 1, parametrizing a continuous family of finger so- 
lutions. However, experiments only showed fingers with 
A = 1/2. This selection problem was understood only 
three decades later by different groups [30l.[3ll.[32l[33l[3^. 
They included surface tension into the boundary condi- 
tion for the pressure p on the interface. This bound- 
ary condition also prevents cusp formation within a finite 
time [m H(| ; this leads to a regularized moving bound- 
ary problem. It was shown, using expansion beyond all 
orders and reduction to a nonlinear eigenvalue problem, 
that in the limit of small surface tension only the fin- 
ger with A = 1/2 is stable. Recently it was found that 
the so-called "kinetic undercooling" boundary condition 
also leads to regularization and dynamical selection of 
the Saffman-Taylor finger with width A = 1/2 for in- 
finitesimally weak regularization 37]. We recently have 
pro pos ed a similar regularization mechanism for stream- 
ers [Hill. 

We therefore have superimposed the Saffman-Taylor 
finger with width L/2 as a solid line on the streamer 
in Fig. [31 The agreement is convincing. In Fig. [51 the 
comparison is further elaborated for three different cases, 
where (b) differs from (a) by the initial density distri- 
bution, and (c) by E^ and L. Here solid lines repre- 
sent stages of evolution of the density model from initial 
transients to uniform translation; they indicate the posi- 
tion y{x) of the maximal charge density for every x. The 
dashed lines are the Saffman-Taylor solution ([5]) with the 
selected width A = 1/2. No adjustments are possible, ex- 
cept for an arbitrary translation of the Saffman-Taylor 
finger along the x-axis. This is chosen to overlap with 
the latest stage of the density evolution which is the at- 
tractor of the dynamics (at a later stage the shape of 
the front stays identical and it moves at constant speed). 
Again the agreement is very convincing. A direct conse- 
quence of this agreement is that we expect the field to be 
enhanced by a factor of 2 immediately ahead of the front, 
and the finger velocity to be 2^^, independently of the 
values of L and Eoo . Indeed we observe that this value of 
the enhanced field is reached when the moving boundary 
approximation is most accurate, i.e. when the width of 
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FIG. 6: Solid lines: contour lines characterizing the simula- 
tions at time steps of At = 100; dashed curves: the uniformly 
translating Saffman-Taylor finger solution Q with A = 1/2. 
(a) background electric field Eoc = 0.5 and width L — 256 
with an initial seed smaller than the steady state solution 
(same as in Fig. [3}, (b) same as in (a) but with an initial seed 
wider than the asymptotic solution (here it takes longer to 
reach uniform translation) and (c) i?oo = 0.6, L = 376. 



the space-charge layer is much smaller than the radius of 
curvature of the front. The former is rather independent 
on L and Eqo, while the latter is of order L, since the 
width of streamer approaches L/2. However, since large 
L also leads to branching, this behavior is observed only 
for parameters slightly to the left of the phase-separation 
curve of Fig. O 



C. Open problems for boundary analysis 

However, this apparently very successful interfacial 
model relies on four approximations. 

1. For the front motion, the electron drift velocity 
v in the local electric field E is increased by a 
diffusion-reaction correction (l9j . The present sim- 
ulations show that the streamer velocity in the 
maximal electric field E + can be linearly interpo- 
lated by v = 1.312 E+ + 6 ■ 1(T 4 within the ex- 
plored field range, giving values closely below the 
analytical result v = \E\ + 2^J D^e" 1 /^ for pla- 
nar fully relaxed fronts Such a velocity cor- 
rection v = c E + can be absorbed completely into 
rescaling time with c. 

2. The streamer interior is not field free immediately 
behind the ionization front as Fig. 2]shows. Conse- 
quently, in contrast to the prediction of the moving 
boundary approximation the front obtained from 
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the minimal streamer model (H)-© is not com- 
pletely equipotential, as Fig. [3] shows. 

3. The space charge layer has a finite width. As a 
consequence, it also can be seen in this figure that 
the electric field is not enhanced by a factor of 2 but 
somewhat less, while the interface position agrees 
very well. 

4. The interfacial approximation breaks down at the 
sides of the streamer finger where the local elec- 
tric field E is too low to sustain substantial ion- 
ization, e _1 /' E ' <C 1, while the interface between 
two fluids in the Saffman- Taylor finger, of course, 
continues along the whole channel length. There- 
fore the mathematical similarity between Saffman- 
Taylor fingers and streamer fingers holds only close 
to their tips, while the analytical construction of 
fingers requires their whole length. 

These observations pose new challenges to the theoret- 
ical understanding of finger selection in moving boundary 
problems. 



properties of a bunch of streamers. We obtained a phase 
diagram spanned by the electric field Eao and the spatial 
period L, see Fig. [5J For L and/or Eoo large enough, 
the streamers branch similarly to single streamers. For 
L and/or small enough, the streamers do not branch 
and approach the width L/2. Furthermore, we used a 
moving boundary approximation to derive surprisingly 
accurate predictions. We showed that close to the brach- 
ing line of the phase diagram, the enhanced field at the 
tip of the streamer is close to 2Eoo, where E^ is the 
background electric field applied between the electrode. 
Moreover we showed that the shape of the front is well 
fitted by the selected Saffman- Taylor finger derived ana- 
lytically from the moving boundary approximation. 

Certainly there are still many open questions about 
this topic. Further investigations should extend our 
model to three spatial dimensions and to a wider va riety 
of media, including nonlocal ionization mechanisms [2l| . 
A rigorous analysis of the problem of finger selection in 
this context of interacting streamers would also prove 
valuable both for the pattern formation community and 
for an improved understanding of streamers. 



IV. SUMMARY AND CONCLUSIONS 

This paper presents, up to our knowledge, the first 
studies on the full dynamics of multiple interacting 
streamers. By using a simplified but physically relevant 
model, we were able to focus on the main effects of the 
interaction and stress the most general electro-dynamic 
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